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Abstract. Motivated by the mean value property of harmonic functions, 
Ps| we introduce the local and global median value properties for continuous 

functions of two variables. We show that the Dirichlet problem associ- 
^ ated with the local median value property is either easy or impossible 

to solve, and we prove that continuous functions with this property are 
1 -harmonic in the viscosity sense. We then close with the following con- 
jecture: a continuous function having the global median value property 
I— I and prescribed boundary values coincides with the function of least gra- 

Ph dient having those same boundary values. 

< 

i 

^ 1. Introduction 

Let Q c be a domain with smooth boundary (9Q, let B(x, r) denote a 
^ ball with radius r, center x, and boundary dB{x, r), and let j^fdji denote the 

O average of / over the set E with respect to the measure yu. It is well-known 

that the continuous function u is harmonic in Q if and only if 

^ C 

QQ u{x) = T u{s)ds whenever x g Q and B{x,r) d Q; (1) 

JdBix,r) 

^ this is the famous mean value property of harmonic functions [4J. As is also 

^ well-known, imposing a Dirichlet boundary condition determines a unique 

• fh solution of the functional equation ([T]). This equation has been well under- 

^ stood for a long time, and it is natural to wonder what happens if we replace 

the averaging operator in ([T]) with a different statistical measure. 

What happens, for example, if we replace the average in ([T]) with the 
median of u over spheres centered at xl Requiring u to be continuous makes 
this question meaningful; as discussed below, medians are not well-defined 
for functions that are merely measurable or integrable. 

When N = 2, which we assume henceforth, we can use elementary tools 
to analyze continuous solutions of the functional equation 

u(x) = median { u(s) } for x e CI and < r < R(x) , (2) 

sedBix,r) 
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where /? : Q ^ R is a continuous function such that 

0<R(x) < dist(x,5Q). 

Since R(x) can be strictly smaller than the distance dist (x, dCl) from x to 
(9Q, equation ([2]) defines the local median value property, a property that 
does not seem to have been studied before. We will show that equation 
([2]) is closely related to a nonlinear partial diff'erential equation (PDE) that 
has a reputation for being difficult. The degeneracy of this particular PDE 
precludes direct applications of standard techniques, making the basic prob- 
lems of existence, uniqueness, and regularity of its solutions rather thorny. 
Our analysis of equation ([2]) clarifies some of these issues: if we impose a 
Dirichlet boundary condition on solutions of Q, we will see that we should 
generally expect either nonexistence or nonuniqueness, depending on the 
interaction of the boundary data and the geometry of (9Q. 

The mean value property ([T]) of harmonic functions is global, in the sense 
that it holds at x g Q for any ball B{x, r) that is strictly contained in Q. The 
global version of ([2]) is also of interest; functions with the global median 
value property satisfy the functional equation 

u{x) = median { u{s) } whenever x g Q and B{x, r) d Q . (3) 

sedB{x,r) 

Combining equation ([3]) with a Dirichlet boundary condition yields a prob- 
lem that is harder to solve than its local counterpart. We conjecture that the 
solutions of this global problem are precisely the functions of least gradi- 
ent, for which existence and uniqueness results are only available when Q 
is strictly convex ifTSll . 



2. The local median value property 

If / is a measurable function with respect to the measure then m is a 
median of / over the measurable set E if and only if 

LliE) LliE) 

li{xeE : f(x) > m] > '-^ and ii{x e E : f(x) < m] > . 

It is easy to see that such a median m need not be unique. If / is integrable, 
then m is a median of / over E if and only if [ITSll 

I - m\ djiix) = min ] I \f{x) - y\ d^(x) 

Je [Je 

but integrability still does not guarantee a unique median. Using this vari- 
ational characterization, however, Noah proved in [10 ] that continuity suf- 
fices to determine median values uniquely. By diff'erent methods, Waksman 
and Wasilewsky proved the same result in [14J for continuous functions on 
planar domains. By working exclusively with continuous functions, we are 
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therefore assured that medians over spheres are well-defined and that the 
functional equations ([2]) and ([3]) make sense. 

To understand ([2]), we first observe that any continuous function whose 
level sets are monotonically ordered straight lines has the local median 
value property. On the other hand, if w is a smooth function such that 
Du(xq) ^ at xo g Q, then the Implicit Function Theorem guarantees 
that, on some sufficiently small ball B(xo, r), the level sets of u are smooth 
curves whose levels vary monotonically; if u also satisfies ([2]), the level 
curve through xq must be a diameter of B(xo, r). The level sets of a smooth 
solution of (|2]) must therefore be straight lines wherever its gradient does 
not vanish. Our first result, Theorem[T} eliminates this smoothness assump- 
tion. Before going through its proof, we note the following obvious and 
useful result, a weak maximum principle for solutions of ([2]). 

Proposition 1. A nonconstant continuous solution o/(|2]) cannot have strict 
interior maxima or minima. 

Figure [T] illustrates a variant of this maximum principle that will be es- 
sential below. It depicts a continuous function u that has a "ridge line," a 
curve along which u = a such that u < aon either side (of course, u could 
just as well be larger than a on either side of this curve). Our proof of The- 
orem [T] relies on the fact that such a situation cannot hold if u has the local 
median value property everywhere. 




Figure 1. A continuous function cannot satisfy (|2]) along a 
ridge line. 

Theorem 1. Let u be a nonconstant continuous solution of (|2]). If K is a 

path-connected component of some level set of u, then dK Pi Q consists of 
line segments with endpoints on dCl. 

Proof Let ^ be a path-connected component of the level set u~^{a), and 
let X e K. The local median value property guarantees that points in K can- 
not be isolated; consequently, at least one nontrivial path must pass through 

X. 
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Step 1. Suppose that, for some r > 0, B(x,r) n u~^(a) consists of 
one simple path F with endpoints y,z ^ dB(x, r). Decrease r (if necessary) 
so that ([2]) holds on dB{x, r). By the continuity of u and the maximum 
principle, y and z must divide dB{x, r) into two open arcs such that u> aon 
one of them and u < aon the other. The local median value property forces 
these two arcs to have the same length; y and z are therefore antipodal points 
and r must be a diameter of B{x, r). 




Figure 2. A simple path through x must be straight. 

Step 2. Now suppose that, for some r > 0, B(x,r) n u~^(a) \ {x] 
consists of simple paths that are not intersected by any other paths. It fol- 
lows from the preceding argument that these paths must be line segments 
emanating from x (Figure |3]). 




Figure 3. Separable paths through x must be straight. 

Continuity and the maximum principle prohibit an odd number of line 
segments from meeting at x, since u must be alternately larger than a and 
smaller than a in the sectors determined by these segments. Moreover, four 
or more segments cannot meet at x : if the angle between two adjacent seg- 
ments were less than or equal to 7r/2, then the local median value property 
would be violated (see Figure]?]). Applying Step 1, we conclude that the 
only path through x inside B(x, r) is a diameter of B{x, r). 

Step 3. If the paths emanating from x can be separated, then the previ- 
ous step shows that, on some sufficiently small neighborhood of x, there can 
be only one straight path through x. In fact, the observations in Step 2 show 
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Figure 4. If 6 were less than then (|2]) would force 
u(y) > a, di contradiction; if 6 were exactly 7r/2, then we 
would have u{y) = a, another contradiction. Consequently, 
e > nil, 

that, if two of the paths emanating from x cannot be separated, then K must 
have a nonempty interior. Consequently, if the interior of K is empty, K 
must be a line segment. By the local median value property, the endpoints 
of this segment cannot be in Q; its endpoints must be on 5Q. The theorem 
therefore holds for path-connected components with no interior. 

Step 4. Finally, suppose that the interior of K is nonempty, pick a point 
y e K, and let z be a point in K such that the segment L connecting y and z 
is inside Q. If L did not belong to K, then the continuity of u would guar- 
antee the existence of a point p on L whose level set component had empty 
interior. By the previous step, that component would be a line segment with 
endpoints on the boundary of Q, separating y from z and contradicting the 
fact that y and z belong to K. This completes the proof. 

□ 

3. The DiRICHLET PROBLEM 

We now seek a continuous solution w of ([2]) that has prescribed values 
along the boundary of Q. Theorem [T] tells us how to go about solving the 
resulting Dirichlet problem: roughly speaking, the level sets of u need to be 
line segments or polygons. 

Figure |5] illustrates a solution of this problem for a nonconvex domain Q 
and very particular boundary data. This example is artificial, however; in 
general, we cannot expect to solve the Dirichlet problem on a region Q that 
is not strictly convex, precisely because we cannot connect arbitrary pairs of 
boundary points with line segments that remain inside Q. See, for instance. 
Figure |6} 

If Q is strictly convex, on the other hand, solving the Dirichlet problem 
is much easier. As an example, let Q c be the open unit ball, and let 
g(x,y) = \y\ for (x,y) e dCl. It is easy to see that, for any a e (0, 1), the 
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Figure 5. A solution of the Dirichlet problem on a noncon- 
vex domain Q. 




Figure 6. The Dirichlet problem cannot necessarily be 
solved if Q is not strictly convex. Here, we cannot connect 
positive boundary values with line segments inside Q. 

function 

( \y\ when \y\>a, 

Ua(x,y) := < Vl - when |_y| < a and \x\ > Vl - , (4) 

[ a otherwise 

has the local median value property and satisfies u = g on dQ. Figure [7] 
provides contour plots of solutions for three values of a, illustrating the fact 
that diff'erent solutions of the Dirichlet problem correspond to the diff'erent 
ways of connecting the boundary values with straight lines. 




Figure 7. Solutions Ua of the Dirichlet problem on the unit 
ball with Ua(x,y) = \y\ on its boundary. From left to right, 
a < V2/2, a = V2/2, and a > V2/2. 
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An existence result for reasonable boundary data follows directly. Specif- 
ically, let Q be strictly convex and bounded, and suppose that g : dQ. ^ R 
is a nonconstant continuous function whose level sets have finitely many 
endpoints. (If g is constant, it follows from the maximum principle that the 
only solution of the Dirichlet problem is constant.) By parametrizing (9Q, 
we can regard g as a continuous 2:7r-periodic function, g : [0, In] R. De- 
fine m and M to be the minimum and maximum of g, respectively, pick any 
A e (m, M), and define A := g~^(A), the set of angles 6 such that g(d) = A. 
Since A has finitely many endpoints, dA = {^i, . . . , Ok] for some k, with 

01 < 02 < . . . < 

Construct a generalized polygon whose vertices on dCl correspond to the 
angles 6i,. . . ,6^. This is a generalized polygon since some of its edges 
could be curved portions of (9Q; this will happen whenever g{6) = Aon the 
closed interval [6j, 0j+i]. As for the other edges, we can proceed as follows. 
Suppose, without loss of generality, that g{6) > A for dj < 6 < dj+i. For 
sufficiently small s > 0, the level set g~^(A -\- s) n [6j,6j+i ] will consist 
of exactly two angles, which we connect with a line segment. If this level 
set never contains more than two angles as s increases, then we are done 
with this sector of Q. Otherwise, there will be a smallest s > such that 
the boundary of g~^(A + s) n[ dj, dj+i ] consists of angles 0i, . . . , 0^ for ^ > 
3. Construct a generalized polygon with vertices corresponding to these £ 
angles, and then repeat this procedure in each resulting sector of Q. The 
compactness of Q guarantees that this process will terminate, yielding a 
solution of the Dirichlet problem; its continuity follows directly from the 
continuity of g. We have thus proven 



Theorem 2. IfQcM? is a strictly convex, bounded open set and g : dCl 

R is a nonconstant continuous function whose level sets have finitely many 
endpoints, then there is at least one continuous function u that satisfies (|2]) 
and equals g on dCl. 



Implementing this procedure for diff'erent values of A can produce diff'er- 
ent solutions; in Figure [7] above, the three solutions correspond to choosing 
A = a for three diff'erent values of a. Consequently, we generally expect to 
find multiple solutions when Q is strictly convex, and an intriguing prob- 
lem for future work is to quantify the number of solutions that correspond 
to a given boundary condition. For a domain Q that is not strictly convex, 
it would also be interesting to characterize those boundary conditions for 
which solutions exist. 
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4. 1 -HARMONIC FUNCTIONS 

For p > \, the continuous function u is /7-harmonic if and only if it is a 
viscosity solution of 

-6iN[\Du\P-^Du) = 0. (5) 

Juutinen et al. proved BUI that this is equivalent to the more common defini- 
tion of /^-harmonic functions as weak solutions of ([5]), the Euler-Lagrange 
equation corresponding to the variational integral 



/ 



\Dv\Pdx for vgWI'^(Q). (6) 



When p > 1, the variational approach to (|5]) is especially elegant, since the 
associated functional ([6]) is strictly convex and weakly lower semicontinu- 
ous and the Sobolev space W^'^(Q) is reflexive. 

Formally setting = 1, one can define a 1 -harmonic function w to be a 
viscosity solution of 

( Du\ ^ 
-div =0, 
I \Du\ I 

an equation that is not nearly as well understood as (|5]). It is more com- 
plicated than ([5]) in several regards; the functional (|6]) is no longer strictly 
convex, for example, and W^'^(Q) is not as hospitable a setting for func- 
tional analysis as W^'^(Q) (as explained further in, e.g., [3J and [5J). 

For various reasons, it is more convenient for us to define to be 1- 
harmonic if and only if it is a viscosity solution of 

-A,„:=-|D„|div(^) = 0. (7) 

As clarified below in Definition [Tj we can loosely interpret this equation to 
mean that either \Du\ = or div (Du/\Du\) = 0, and we will see that contin- 
uous solutions of ([2]) are 1 -harmonic. We first define viscosity solutions of 
(|7]), referring to [[Sjl for more details when p > I. 

Definition 1. The continuous function S : Q ^ R is a viscosity superso- 
lution of (|7]) if and only if 

(1) u f oo and 

(2) whenever xq e Q and (p e C^(Q) satisfies 
u(xo) = (p(xo) , 

u(x) > (p(x) for X ^ xo, and (8) 
D^(xo)^0, 

ip also satisfies 

-Ai^(xo) > 0. 
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The continuous function w is a viscosity subsolution of (|7]) if and only if -u 
is a viscosity supersolution of (|7]), and u is 1 -harmonic in the viscosity sense 
if u is both a viscosity subsolution and a viscosity supersolution of (|7]). 

The work in [8J motivates this definition, as Juutinen et al. proved that, 
when p > 1, we only need to work with test functions whose gradients do 
not vanish at the point xq of interest. It then follows that our definition of 
1 -harmonic and the usual one are equivalent. 

Theorem 3. A continuous function w : Q ^ R that satisfies Ij^ is l- 
harmonic in the viscosity sense. 

Proofi Let xq g Q, and let ip g C^(Q) be a test function that touches u from 
below at xq. Using ([2]), ([8]) and the monotonicity of the median, we have 



so that 



median { (f){s) ] < median { u{s) } = w(xo) = v^(xo) . 

dseBixo,r) sedBixo,r) 



(p(xo) - median { (p(s) } > . (9) 

sedB(xo,r) 



Identity (11) in [6 J establishes that 



(p(xo) - median { (p(s) ] = -— Ai^(xo) -\-o(r), 

sGdB(xo,r) 2 

SO dividing by and sending r ^ in ([9]) shows that u is a viscosity 
supersolution of (|7]). Similarly, uis a viscosity subsolution of (|7]), proving 
that u is 1 -harmonic in the viscosity sense. □ 

The key result above, identity (11) of [O, is the reason for including the 
pref actor \Du\ in our definition of the 1-Laplacian Ai. Not coincidentally, 
our definition of Ai coincides with the elliptic part of the operator in the 
widely studied mean curvature equation, 

/ Du \ 
.,-|fl.|div(— 1 = 0. 

Previous work relating this parabolic equation, medians, and median-like 
operators (Q, 113, 13, 0,113) motivated the present work. 

5. Functions of least gradient 

In a series of papers summarized in ifTSll . Ziemer et al. analyzed func- 
tions of least gradient, i.e., functions v e BV(Q.) Pi C(Q) with minimal total 
variation and prescribed values along dCl. They proved, in particular, that a 
unique function of least gradient exists if Q is strictly convex and the bound- 
ary function g : ^ R is continuous. One can assemble this function of 
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least gradient by constructing its level sets: they are minimal surfaces whose 
boundaries are determined by the levels of g. 

This, of course, is exactly how we constructed solutions of the Dirichlet 
problem earlier, since one-dimensional minimal surfaces in the plane are 
lines. Our method did not produce a unique solution, however, as we only 
considered local conditions; the global constraint of minimizing total vari- 
ation renders functions of least gradient unique. It seems reasonable, then, 
that there should be a strong connection between the global median value 
property ([3]) and functions of least gradient. To that end, we propose 

Conjecture 1. Suppose that Q is strictly convex and that g : (9Q R is 
continuous, and let u* be the function of least gradient on Q that equals g 
on (9Q. 

(1) There exists a unique continuous solution uoflj^ such that u = g 
on dQ. 

(2) u = u\ 

It is easy to see that this conjecture holds in certain cases. For example, 
among the family {Uq.} of solutions of ([2]) defined by equation (|4]) earlier, 
only ^1/^2 has the global median value property. It follows easily from the 
coarea formula that this function's total variation is smaller than that of any 
other Ua, and it is indeed the function of least gradient on the unit circle 
with this boundary data. 

Progress on this conjecture (and on other questions related to the results 
of this paper) will likely benefit from Juutinen's work [TJ on /^-harmonic 
functions and functions of least gradient. 
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